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Abstract 



< 

. We study the mathematical properties of a general model of cell division structured with several 

' internal variables. We begin with a simpler and specific model with two variables, we solve the 

eigenvalue problem with strong or weak assumptions, and deduce from it the long-time conver- 
gence. The main difficulty comes from natural degeneracy of birth terms that we overcome with a 
regularization technique. We then extend the results to the case with several parameters and recall 
the link between this simplified model and the one presented in ; an application to the non-linear 
problem is also given, leading to robust subpolynomial growth of the total population. 
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'jh ■ 1 Introduction 

1.1 Presentation of the Model and Link with Other Models 

Regulation of the cell division cycle governs the development of all organisms. Understanding it is 
central to the study of homeostasis, tumour growth and cancer, but is made particularly difficult due 
to the numerous phenomena that can have an influence on it (see for instance and [13] , or [3D] for 
a general presentation of the cell cycle). 

For these reasons, several models have been proposed usually structured by a single variable (age, 
size, etc) [16]. But modern biology offers more accurate structuring variables as proteins or molecular 
content [35] . 

In order to investigate qualitatively the long-time behaviour of a cell population, we consider here 
a general model structured both in age, represented by the variable a, and in another aggregated 
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variable x. Namely: 



-Qin+ -^\r{a,x)n\ + B{a,x)n = 0, a ^ 0, x ^ 0, 

n(t, a = 0, x) = 2 J b{a, x, y)n{t, a, y)dy da, (1) 

r+(a, X = 0)n(t, a, X = 0) = V a ^ 0. 

This variable x can represent one of the various proteins produced (cf. for instance [36] or [Q), the 
maturity of the cell (as in [1] and [2] for instance), its size, its DNA content (as in [3] for instance) etc. 
Our study can be generalized to the case when x € , that is, when several phenomena influencing 
the cell cycle are taken into account (see part I4.ip ; in part 14.31 we also investigate the possible 
application of this model to a non linear two cell-compartment model (to model proliferating and 
quiescent cells). 

We have supposed here that age evolves like time, i.e. ^ = 1- The function T = ^ represents the 
rate at which the x content of a cell increases with age. The function B represents the total division 
rate, and 6(a, x, y) is the repartition function of a mother cell of y— content to a daugther cell of x— 
content. We impose for consistancy (see [B]) 

Bia,y) = J b{a,x,y)dx, (2) 

yB{a,y) = 2 J xb{a,x,y)dx, b{a,x,y) = b{a,y - x,y), b{a,x > y,y) = 0. (3) 

It means conservation of the number of cells, conservation in the x variable and symmetry of the 
division, when x represents a molecular content. 

We have taken the coefficients independent of time (see |14) for a model with time-periodic coeffi- 
cients). If we suppose that the coefficients T and B do not depend on x and that for the solution Tn 
vanishes for x = and x = oo, integrating equation ([1]) in x and denoting J\f{t,a) = J n{t,a,x)dx 
gives: 

'' §iJ^+-§-J^ + B{a)M = 0, a^O, 

(4) 

M{t,a = 0) = 2/ B{a)N{t,a)da. 

This is the classical linear McKendrick-Von Foerster equation, with a death term which is exactly half 
the birth term (see [31] for a complete study of this equation in the case whithout death term, and 
|25j for extension of this equation to a non- linear case) . 
If we suppose that the coefficients F, B and b do not depend on the age variable a, that the integral 

oo 

N{t, x) = J n{t, a, x)da converges and that lim n(t, a, x) = 0, we can integrate ([T]) in a and find: 







iM+£ [T{x)M] + Bix) M{x) = 2 / 6(x, y)M{t, y)dy, 

(5) 

F(x = 0)7V(t, X = 0) = V t ^ 0. 



If F = 1, we find the pure size-structured model, which has been studied in [M], [2H], [32], [M] for 
instance. In the case F = 1 or F = x'^, existence of a solution to the eigenvalue problem for a general 
b is proved in [23] using approximation scheme. 
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A model describing the dynamics of a cell population divided into proliferative and quiescent com- 
partments was presented in [6j and [7j . It can be written: 



P+^P+ x)p\ + {B{a, x)+di + L{a, x)) p - G{N{t))q = 0, a ^ 0, x ^ 0, 



at 

§iq + (G[N{t)) + d2)q = L{a, x)p, a ^ 0, a; ^ 0, 

(6) 

p{t, a = 0, = 2 / b{a, x, y)p{t, a, y)dy da, 

p, q ^ 0, J p{t, a, x) + q{t, a, x) dx da = N{t). 
We will see in part 14.31 how it can be reduced to the study of problem ([TJ . 

1.2 The Eigenvalue Problem 

In order to study the asymptotic behaviour of the solution of problem ([1]) , we consider the eigenvalue 
problem (see [6j and find {\q,N) solution to 

+ £[T{a,x)N] + {Xo + B{a,x)) N = 0, a^O, x^O, 
N{a = 0, x) = 2 / b{a, x, y)N{a, y)dy da, (7) 

iV ^ 0, jNdxda = l. 

This is an original problem which can be seen as a usual Cauchy problem where the initial data is 
related to the "future" . 
It is useful also to study the adjoint problem: 

-fj - r(a, x)-§^^ + (Ao + B{a, x)) = 2 / 6(a, y, x)(/>(0, y)dy, a ^ 0, x ^ 0, 

(8) 

'^0, j(l)Ndxda = l. 
Here we make the assumptions: 

r(a, x) = xr(a, x), (9) 

with _ 

r(a, x) > for X « 0, 

(10) 

r(a,x) ^0 for X ^ xjv/- 

Then there is no need of boundary condition at x = and conservation in the x variable is enforced 
according to the biophysical interpretation, when x represents a molecular content. The fact that 
r becomes negative beyond a maximum value xm means that the x— content of the cells remains 
bounded, but the results can be generalised to the case when T remains nonnegative everywhere. 
As a consequence of condition integrating equation ([7]), we have for vanishing at infinity 

Ao = / B{a,x)N{a,x)dx da, 



3 



(in words, the population number can only grow by cell division). Integrating the equation ([7]) against 
the weight x, as soon as lim f xN{a, x)dx = 0, we have, using ([3]) 

a— >oo 

XolxNia,x)dxda = lTia,x)Nia,x)dxda, 

(in words, the total molecular content can only increase by the reaction term T). Integrating the 
equation against the weigth a, as soon as lim J aN{a, x)dx = 0, we have 



A, 



t j aN{a,x)dxda + J aB{a,x)N{a,x)dxda = 1. (11) 



Also, we can do as in [24] and for < r/ < 1, if lim f e'^'^^"'N(a,x)dx = 0, multiplying by e'*'"^" and 

a— »oo 

integrating, we find 

V < < 1, / e^o'^'^Nia, x)dxda ^ ^ 



1 — r] 



We can reduce the study to the solutions on the domain (a, x) € M+ x [0,XAf]. Indeed, using the 
method of characteristics based on the solution to the differential system parametrized by the Cauchy 
data (a, x) : 

£X{a,x) =T{a,X{a,x)), a^O,x^O, 

(12) 

X{0,x)=x, x^O, 

Cauchy-Lipschitz theorem gives us, as soon as F E C^(M_|_ x M_|_) for instance, the existence and 
uniqueness of the flow X(a,x). We denote Y{a,x) the inverse flow, defined by 

Y{a,X{a,y)) =y. (13) 

We deduce from ([10]) that for all x ^ xmi for all a ^ 0, X{a,x) ^ X(a, xm) ^ xm- The formula 
(|19p . proved below in lemma [LTl shows that if the solution N verifies N{a = 0, x ^ xm) = 0, then 
A^(a ^ 0, X ^ X(a,a;A/)) = 0, and in particular N{a,x ^ xm) = 0. Thus, we add the following 
condition to problem ([7]): 

iV(0,x^XAf) = 0, (14) 

and it allows us to restrict our study to the compact set [0, xm]- We could also exchange assumption 
([9]) with the following one: 

N{a ^ 0, X = 0) = 0, b{a ^ 0, 0, y) = 0. (15) 

Contrarily to the solution A^, the solution of the adjoint problem ([8]) does not necessarily have its 
support in [0,xm]) but we make the following assumption on the support of the function b (cf. the 
proof in section I2.3P : 

b{a, X, y ^ X{a, xm)) = 0, (16) 

which implies that 6(a, x > XM,y) = 0, and thus, by formula ([T7]) proved below, it implies ([H]) . 

In all the following, if nothing is specified, we suppose assumptions ([9]), ([TO]) . ([2]) and ([3]) are satisfied. 
We denote X{a, x) the characteristic flow solution of (I12p and Y{a, x) the inverse flow defined by (I13p . 
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1.3 Reformulation of the Problem with the Method of Characteristics 

We first give the following formulae, on which the proofs are based. 

Lemma 1.1 ForT £ C^(M+x[0, j;a/])) under assumption a solution N to ^ verifies the following 
formula (as soon as the integral converges): 

N{0,x)=2 b{a,x,X{a,y))N{0,y)e <' dyda, (17) 



and we can write: 

- f I >'0 + 4-'^(s,X(s,x))+b(s,X(s,x)) \ds 

N{a,X{a,x)) = N{0,x)e /, ^ ° ^ ' ^ ' K' »S ^ ^^^^ 

or also 

-f{Xo+f^r{s,x{s,Y(a,y)))+B{s,x{s,Yia,y)))}ds 

N{a,y) = N{0,Y{a,y))e « . (19) 

Under assumption U6\) a solution (p to ^ verifies (as soon as the integral converges): 

f f , , -J {Xo+B(s,X{s,x))}ds 

(P{0,x)=2 J J b{a,y,X{a,x))^{0,y)e ' dyda, (20) 



and as soon as the integral converges, we can write: 

OO Xm s , , 

r r , , -\Q{s-a)-f B[a,X{ay{a,x)))du 

cPia,x) = 2j J b{s,y,X{s,Y{a,x)))my)e ^ ^ ' dyds. (21) 

a 

~ / {Ao + ^r(s,X(s,x))+B(s,X(s,x))}ds 

Proof. We set N{a,x) = N(^a,X{a,x))eo ^ and rewrite equation ([7]): 

3 ~ 8 8 f \ Xo+-^r(s,X(s,x))+B(s,X(s,x)]ds) 

lN{a,x)=[lN + l[TN] + {Xo + B)N]{a,Xia,x))J^ a. I ( I ( U 

This gives the equalitities (jlSp and (|19p . We can rewrite the boundary condition of problem d?]) as: 

y - 1 {\o + ^T{s,x{sX{a,y)))+B(s,x{sX(a,y)))\ds 

N{0,x) = 2 bia,x,y)N{0,Y{a,y))e i ^ ' ^ ' ^ ' ''""'^ ^' ^' ^'"^^^^ dyda, 





-/^r(s,X{s,y))ds 

or equivalently, changing variables y — > Y[a,y) and noting that e ° = ^-^(ay)' 

formula p!7|) . 
In the same way, we set 

~ - J {Xo+B(^s,X(s,x))]ds 

cl){a,x) = cl)[a, X{a,x))e ° 
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and rewrite equation ([8]) as 



a 



-( f \\o+B(s,X{s,x))}ds) 

m^Ka,x) =[9-^ + T^cj>-{Xo + B)4>]{a,Xia,x))e 

-} {X()+B{s,X(s,x))\ds 

= -2 J b{a,y,Xia,x))cPiO,y)e ' ' dy. 



Integrating along a on M^. we get formula (j20p , and integrating from to a we then find formula ([21 



1.4 Some Examples 

Case 1: for the function F, we can take for instance T{a,x) = Cix(xm — x)- We easily calculate 
that 



X(a, x) 



X ■ XMG' 



Xm + xie'^MCia _ X) ' 

so for X ~ 0, one has X{a,x) ~ xe^'^"'-^^"'. 

Case 2: another possible example is to take T{a, x) = Cix'^{xm — x)^, with < a < 1 and /3 > 0. 
For X > 0, we then have 

Xia,x) «,«o (Ca + x^-")T^. 
Case 3: biological considerations lead [5j and ^ to take for T (case illustrated in fig. [1]): 

X Co 

F(a,x) = ci- (ri - r2e~''^'') - C2X, n, r2 > 0, — < n - r2. (22) 

1 + X ^ ' Cl 

In this case, xm = fj-n — 1 and F(0, xq) = with xq = ^{"i^i — ^2) — 1 > 0. 




Figure 1: Form of F defined by (|22p. and its related characteristics. 



For the division rate B, we can take 

B{a,x) = C2X^lA*^a^A^, 7^1, O^A*<Ai^oo. 



(23) 
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To define the repartition function b (see [22], [37] and [38] for biological motivations) we can choose a 
uniform repartition: 

b{a,x,y) = l^^y . 

On the opposite, the classical example of equal mitosis is given by 

b{a,x,y) = S^^RB{a,y). 

It is the case for instance if x represents the DNA-content of the cells. In this last case, formulae p7|) 
and (pOj) of lemma 11.11 are no longer valid; see part 13.31 for an adaptation of the proof. A study of 
the link between the repartition function b and the proliferation rate Aq, for the pure size-structured 
model, can be found in |26j . 

We can also generalize the uniform repartition by b{a,x,y) — jj(jT:^^jj^r]ys^x^{i-ri)y with < r/ < ^, or 
take a Poisson distribution for b. 



2 Resolution of the Problem in a Regular Case 
2.1 Main Results 

In this part, we consider regular functions b and B in order to obtain strong compactness and regular 
solutions with non extinction, i.e., Xq ^ 0. We need the following assumptions. 

rGC^(M+,C2(M+)). (24) 
This condition is used to prove the regularity of the solution N (else it can be weakened: see part 

15:21). 

B eC{R+x[0,XM]), ^B €C{R+x[0,xm]), B{a,0) = 0, (25) 

beL^iRl), b{a,.,y) £CI{{x(£[0,xm]; x^y}). (26) 

The function b is not necessarily continuous for x G M+ since b{a, x > y,y) = 0, but it is necessary 
to suppose it is regular for x ^ y in order to prove strong compactness. On the previous examples, 
except the case of equal mitosis, and with 7^1, the function b verifies assumption (j26p . 



oo xm " / \ 

- J B[s,X(s,x))ds 

e dx da <oo. (27) 



This is a key assumption, which is used to obtain compactness and a solution N in L-^(IR^). We can 
give an equivalent formulation in terms of partial differential equations, as follows. 
Let V the solution of the following problem: 

+ ^ [r(a, x)v] + B{a,x)v = 0, 0, xm > x ^ 0, 



v{a = 0, x) = l O^x^ Xm, 



(28) 



then ()27p means v G L^(IR^). Assumption (j27p can be replaced by slightly more general ones, allowing 
in particular J B(^a, X{a, x)^da to be finite: see for instance paragraph! 
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Finally, we need the following two assumptions to prove uniqueness of a solution. 

3^1 ^ ^* > 0, V^* ^ a ^ yli, Vy G ]0,xm[, B{a,y) > 0, (29) 

VxG]0, XAf[, Vy G]0,XAf[, j b{a, x, X{a,y))da > 0. (30) 
Assumption ([30]) implies that Suppa{b) is unbounded and that lim X(a,0 < y < xm) = xm- 

a— >oo 

Uniqueness could also be proved under other assumptions than (p9|) and ([30|) : see for instance part 

EH 



In the previous examples, the function F verifies assumption (j24p but we need to regularize B in order 
to obtain (I25p . and also to regularize b, in the case of equal mitosis, to ensure (j26l) . The positivity 
conditions (j29p and (jSOp stand if = +oo in the definition of i?, but assumption ()30p is not true in 
the case of equal mitosis. It remains to check assumption ()27p . If F(a,x) = Cix{xm — x), we have to 

prove the convergence of J J e o dxda, or equivalently of 

J J J e'^''^" dxda = j j dxdu = j j ^-j-^dXdu. 



Since 7 > 0, this integral converges, which gives us assumption ([2] 
If F(a,x) ~xKiO Cx"' with < a < 1, we have to study the convergence of 



00 



1= e dxda^C e-" da. 



1+ 





Since 1 + > 0, this integral converges. 

Theorem 2.1 Under assumptions p4\ ), [24\)-(27\), there exists a solution N € C(J(M^) to ^ for a 
Xq > 0. Moreover, N(a,x ^ X{a,XM)) = 0. 

Theorem 2.2 Under the assumptions of theorem \2.1\ with the additional assumption there 
exists a solution € C^(M^) of ([^j for a Aq > 0. 

Theorem 2.3 Under the assumptions of theorem \2.1\ and with the additional assumptions I j6]] . ^29\) 
and [3D\). the solutions (Aq, N, cj)) of the eigenproblems ^ and ^ are unique. 



2.2 Regularised Problem and Method of Characteristics 

In order to apply Krein-Rutman theorem (refer to [15J for instance), we first have to consider a 
regularised problem, where the operator is strictly positive. We write 

be{a,x,y) =b{a,x,y) + —, Bs{a,y) = B{a,y) + e, (31) 

XM 



8 



for e > 0. The functions and Bf, verify the same relation ([2]) than b and B. We consider the 
regularised problem: 

^ i^N, + £[T{a,x)N,] + {X + Be{a,x))N, = 0, a ^ 0, ^ x < xm, 

Nir{a = 0,x) = 2 J b^{a,x,y)N^{a,y)dy da, O^x^xm, (32) 

Ne^O, fN,dxda = l, 
and its adjoint: 

-S^^e -r(a,j;)^(/)e + {X + Be{a,x)) (t)e = 2 f cl)e{0,y)be{a,y,x)dy, ^ a, ^ x ^ xm, 
^ 0, J (psNsdxda = 1. 

(33) 

We can write the formula ()17p for this regularised problem: 



, , (-/{A+_Be(s,X(s,j/))|ds) 

N,{0,x)=2 be{a,x,X{a,y))N,{0,y)e' ' 'dyda. (34) 





To find a solution to this equation, we consider the Banach space 

E = cI{[o,xm],R), \\f\\E ■■= ll/IU- + 

and study the linear integral operator defined on E by: 

r -J \\+Bjs,X{s,y))]ds 

g{if)ix)=2jb,{a,x,Xia,y))fiy)eo' dyda. (35) 

Lemma 2.4 Under the assumptions of theorem \2.1\ for all e ^ 0, A ^ 0, the operator is compact. 

Proof. In order to apply the Ascoli theorem, we evaluate, for an arbitrary f E, ^ 1, the 

quantities A(xi,X2) = t/^(/)(xi) - C/f (/)(x2) and A'(xi,X2) = ^^a(/)(^i) " S^a(/)(^2), and prove 
they tend to zero uniformally on xi, X2 when \xi — X2I tends to zero. Supposing xi < X2, we can write: 

00 a . , xm yia.,X2) 

r ~ f B[s,X(s,y))ds ( r f } 

A<2/e < \b{a,xi,X{a,y)) -b{a,X2,X{a,y))\+ / b[a,xi,X{a,y)) \dyda. 

Y{a,X2) Y{a,xi) 

For each term, thanks to assumption (j27p and to the fact that b G L°°(M5|_), we first restrict the 

00 

integral to a ^ A where A is chosen such that the rest J be sufficiently small independantly of xi and 

A 

X2. The first term may be bounded by C\xi — X2I thanks to (j26p . The second term may be bounded 
as 0(|xi — X2I) since Y{a,x) is equicontinuous on [0,A] x [0, xm]- We evaluate A'(xi,X2) in the same 
way, with the help of Ht^/H ^ /, ,\. This proves the equicontinuity of the family gUf), for f G E 

"ax "L°°[[0,xm\) 

and 1 1/1 Is ^ 1, therefore Qf^ is compact. 
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Lemma 2.5 For all e > and A ^ 0, under the assumptions of theorem \2.1\ there exists a unique 
/^A,e > and a unique N^^ e E, N^^> 0, such that 

gUNl) = f,,,,Nl, \\nI\\e = i. 

Moreover, A — > ^A,e is a continuous decreasing function, vanishing when A tends to infinity and taking 
the value 2 when A = 0. 



Proof. The compact operator is strictly positive on E, thanks to the fact that e > 0, so Krein- 
Rutman theorem gives existence and uniqueness of /iA,£ > and A'^^^ > 0. The function A — > ^a,£ 
decreases because Qf^ decreases with A. It is continuous thanks to the uniqueness of the eigenvector 
and to the compactness of the family of operators Qy For f E, f > and ||/||£; = 1, we have: 

gi{f){x) ^ 2 1 f{y)e-''^\M^dyda = ^IHU + ^^/xm J ^^^^^^^ 
We evaluate this quantity in / = A'^^ ^ , and integrate in x : 



2 6 ooXM + 2£ 

/^A,e ^ T , (36) 



A 



so lim /iA.e = 0. Integrating in x, we find 

A^oo 



r r . d / ~f Bjs,X(s,y))ds\ 

J gi{f){x)dx = 2j f{y)e-^--^-e jdyda. 
Integrating by parts, we get: 

r r ,^ . - f B(s,X{s,y))ds 

gf(/)(x)dx = 2 / /(y)(iy-2 / /(y)Ae-(^+^)'^e dyda, 



so for A = 0, / g\{f){x)dx = 2 J f{x)dx, so /io,e = 2. 



2.3 Proofs of the Main Theorems 

Thanks to lemma [231 foi^ e > 0, we can define A^ > as being the only A > such that fJ-x^,e = 1- 
We denote A''^ the associated eigenvector, with IIA'^PH^; = 1. Since = QxA-^e)^ family (A'^P)o^£^i 
is compact in E. Indeed, as in lemma [2^ we can show that it is an equicontinuous family of functions 
and apply the Ascoli theorem. Thanks to inequality (I36p . we know that A^ is bounded so we can 
extract a subsequence iXe,N^) tending to (Aq, A^°) € M+ x A^o ^ 0, Aq ^ and \\N'^\\e = 1- We 
have 

r F -1 {\^+e+B{s,X(s,y))\ds 

N^{x) = 2 / {b{a,x,X{a,y))N^{y) + —N^{y)}e " ^ dyda. 
Under assumption ()27p . the second term in this expression vanishes when e ^ 0, so we have: 

r - l \ >'0+B(s,X(s,y))]ds 

N'>{x) = g'i^{N'){x)=2jb{a,x,X{a,y))N\y)eo'- dyda. (37) 
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We define N by 

n ^ -jUo + ^T'(s,x(s,Y{a,y)))+B(s,x(s,Y{a,y)))}ds 

Nia,y) = aN\Y{a,y))e 0^ V V (38) 

To state tliat tlie so-defined function N satisfies d?]), it remains only to prove that J Ndadx < oo, then 

-J {Xo+B[s,X{s,y))}ds 

we choose a such as / Ndadx = 1. We have / N{a,y)dady = a f N^{y)e " dyda. 
Thanks to assumption (I27p and since N^ £ -L°°(]R_|_), this integral converges. As seen in part 11.21 

integration of equation ([7]) gives Aq = / BNdxda > 0, because if it were zero, that would imply 

xm 

J N°{x)dx = 2/ BNdxda = 0, so iV° = 0, which is absurd since HA^'^H^; = 1. The continuity and 



the continuous derivability of N at the points y = X(a, xm) come straigthforward from formulae (|37p 
and (jSSp . and from the assumptions of regularity ()24p ~ (j26p : it ends the proof of theorem I2.1[ 

In the same way than for the direct problem, formula (pO|) leads us to study the integral operator 
defined on the Banach space E by 

r r -( f i Xo+b(s,x(s,x)) \ds) 

Gll{f){x) = 2j J b{a,y,X{a,x))f{y)e 'o^- ' dady. 



The previous study for can be carried out in the same way, and we find an eigenvalue Ai > 
and an associated eigenvector (fp ^ E for the operator Q^*- We choose H^^He = 1, and define (/)(a, x) 
by the formula (|2ip . It satisfies ([8]), and it only remains to check that J" (pNdadx < oo. We have: 

oo Xm oo Xm 

J (PNdadx = 4a j j j j (l)^{y)b{s,y,X{s,Y{a,x)))e-^'^-^''^'-''^ 

a 

- f B ((T,X (a,Y(a,x)) ) da- f da:r(a,X (a,Y(a,x)) ) da 

N'{Y{a,x))e i ^' ^' i K' K' ^' ») ^xdsdyda. 

Changing variables x — > Y{a, x) as we did in part 12.21 and integrating according to the s— variable 
before the a— variable, one gets: 



oo Xm Xm 



/,iV.<,.. = 4„////(.)iV»,.)K.,.,X,...))( 



g-Aos _ g-Aig -/B(a,X(fT,a;))da 

■)e dxdyds. 



Ai — Ao 



This integral converges since Aq, Ai > (we do not know yet if Aq = Ai; if it is the case we just replace 
- — ^^Z^o se""**"*). We choose a such that J (pNdxda = 1. It ends the proof of theorem 12.21 

For uniqueness, assumption (f30]l implies that for all x G ]0,a;M[, we have N{0,x) > and </>(0, x) > 
0. Thanks to Fredholm alternative, uniqueness of a solution to problem ([Tj) implies that to problem 
dH]), so we only have to prove uniqueness for one of that two problems. Consider two solutions (Ai, A'^i) 
and (Ao, (p) of ([7]) and ([8|). Then we have: 



-Ai / Ni(pdxda = J y-§^Ni + ^[r{a,x)Ni] + B{a,x) Nij (pdxda 

-- j(^--§^^-T{a,x)^4> + B{a,x)^-2 J 4>{0,y)b{a,y,x)dy^ Nidx da = -Ao / Ni(l)dxda. 
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Thanks to assumption ([30|) . we have / Ni(j)dxda > 0, and it impHes Ai = Aq. 

Let us now suppose two solutions Nq and Ni of problem ([7]) for the same eigenvalue Aq- We use the 
following lemma (for a proof, we refer to [3T], proposition 6.3.) 

Lemma 2.6 Let f G 5'(IR^) be a solution of 

^/ + ^[r(a,x)/] + (Ao + 5(a,x))/ = 0, o ^ 0, a; ^ 0, (39) 
then I/I satisfies the same equation. 

We set = — A^i-I According to lemma [2^ verifies equation ([39]) . We take ^ as a test function: 

(j) {Q,x)b{a,x,y)\N — Ni\{a,y)dx dy da = 2 J (j){0,x)\ jj b{a, x,y){N — Ni){a,y)dx dy \da. 

Thanks to the fact that </>(0, x) > on ]0, xjv/i, it implies that b{a, x,y){N — Ni)(a,y) is of constant 
sign. We integrate in x and find that B{a,y)(N — Ni){a,y) is of constant sign. Using assumption 
(j29|) . — A^i is of constant sign on [j4*, yli] x]0, Formula of characteristics (fT9]) then gives that 
{N - Ni){0,]Y{A*,0),Y{A*,xm)[) is of constant sign. Since Y{A*,0) ^ and Y{A*,xm) ^ xm, we 
deduce that it is of constant sign on {0}x]0, xm[, thus, using once more the formula of characteristics, 
on R\. Since Jf{N- Ni)dx da = 0, we deduce = A^i. 

3 Various Generalizations of the Results 
3.1 Division Rate With Compact Support 

We have used previously strong assumptions on the support of b and B. Here we relax it and suppose: 

3A>0, Suppaib, B) C [0, A]. (40) 
To illustrate this on a simple example, let us take, for i? > constant, B{a, x) = B\a^A- We integrate 

XM 

formula (|35p and find, for any f E such that J f{x)dx = 1 : 





XM A XM 

J gx{f){x)dx = 'iB j j /(x)e-""-^'^dxda = ^^(1 - e-("+^)^) = /.a- 





It implies that there exists A > such that = \ iff fi\=o > 1, that is iff B > ^Log{2). 
As seen on this example, we need a condition to ensure that Aq ^ 0. Namely: 

-fB{s,X{s,x))ds 

j b{a,x,y)e Q dx 1 

a = sup — ^ < -. (41) 

a,y B{a,y) 2 

(We can divide by B{a, y) without loss of generality since B{a, y) = implies b{a, x, y) = 0.) 

Such a condition is not surprising and occurs always in structured population equations (see [31]). 
We check easily that it generalizes the condition of the example B{a, x) = Bla^A- 
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To have uniqueness of a solution and of a A > 0, we need to replace assumption (pOj) . Indeed, we have 
noticed in part 12.11 that it generally implies Suppaih) unbounded. If we suppose X{A^xm) < xm-, for 
XiA,XM) < Xm, if y < mill X(a,XM), then for all ^ a ^ yl, one has X(a,y) < X'^(A,xm) ^ 

X(A,xm) < X, so pO]l cannot be verified. We make the following restrictive assumption: 

VO<a<XM, VO<x<a;M, r(a,x)>0. (42) 

Without assumption (|42|) . we are unable to prove uniqueness of a solution. Refer to the proof of 
theorem 12.31 the fact that J f Ncj) dxda > plays a central role, and here we would only be able to 
prove that > on an interval ]0,xl[ and (/> > on ]xl,xm[ where xi is such that X{A,xl) = xl 
(this proof can be done recursively on the basis of formulae (jl7p and (|20p). We complete assumption 
(f29l) by an assumption on Suppa{b) : 

y x e]0,XM[, y y e [x,xm[, b{A* ^ a Ai,x,y)da > 0. (43) 

Theorem 3.1 Under assumptions ffgj) . (24^- [M\) . (J^ , (JJ^ , there exists solutions (iV, 0) G C^(R^) 
respectively to problems ^ (B^ for a Xq > 0; and Suppa4> C [0,^]. Under the additional assumptions 
[2^) . (4W o-f^d l4^ , the solution (Aq, N, 0) of the eigenproblem is unique and Aq > 0. 

Proof. We proceed in the same way than in part [21 so we let the details of the proof to the reader. 
We first take e > and consider the following regularization of the parameters, which is the symmetric 
of definition ([3T]) : 6e(a, x, y) = b{a, x, y) + ^l^^^, Bs{a, y) = B{a, y) + el.a4.A- With this definition, 
we carry out the same calculations and define the integral operator by the same formula (j35p . 

Lemma 3.2 Under assumptions [W\). [24\)- [26\) . (JO^ j for a// e ^ 0, A € M, the operator is 
compact. 

The proof is the same than for lemma 12. 4[ even simpler since we do not have any difficulty with the 
convergence in the a— variable of the integrals. For this reason, we can now let A < 0. 

Lemma 3.3 For all e > and A G M, under assumptions of lemma\3^ there exists a unique /i^.e > 
and a unique N^^^ G E, N^^^ > 0, such that ^A(^A,e) = MA.eA^A.e' II^A.Jb = 1- ^^^^^ ^« ^ore, 
A fix,e 0' continuous decreasing function which vanishes when A +oo. For A = 0, under 
assumption ( (^i[ ), we have /io,e > 1 for e small enough. 

Proof. The proof is the same than for lemma [231 except that we can let A < 0. It only remains to 
prove that //o,e > 1- We denote by A^q^^ G E the (unique up to a multiplicative constant) solution of 
Gq{Nq ^) = ^o,£-^oe- Defining Nq^s by the characteristic formula ([TOll with Nq ^, and be instead of 
N, B and 6, the function A'^o.e is then solution to the following problem: 

' £No,e + £[T{a,x)No^s] + Be{a,x)No,e = 0, a^O,x^O, 
' l^o,eNo,e{a = 0,x) = 2fbs{a,x,y)No^e{a,y)dyda, (44) 

^ No^e>0, f No,e{a = 0,x)dx = 1. 

Notice that we need e > in order to have convergence of the integral f No^^dadx : indeed, SuppaNo^^ 
is not compact. We integrate this equation in a between and A and in x, and find: 

/ Nos{Ax)dx + {1 —) I Be{a,x)Nos{a,x)dxda = 0. 

J ' /^0,e J 
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Since / B^No^edxda = it can be written / No^e{A,x)dx = 1 - 

Hence, /xo,e > 1 iff f NQ^i;{A,x)dx < i. The characteristic formula allows us to write: 



/ 



NQ^s{A,x)dx= No,e{0,Y{a,x))e » dx 



- J B!r(^s,X{s,x))ds 2 f - J B!r(^s,X{s,x))ds 

No^s{0,x)e dx = be{a,x,y)No^s{a,y)e ° dxdyda. 



Under assumption (I4ip . which is also verified by bs and Bs for e small enough (to be more precise, if 
e~^^ ^ |), it implies / No^i;{A, x)dx ^ j^k^^ ~ I' ends the proof of lemma [331 

We are now ready to prove theorem 13.11 Lemma 13.31 gives us a unique solution and a unique 
Xe>0 for which fix,,e = 1 and ||A^°||£; = 1. 

The family {N^, Ae)o<£<i is compact, so we can extract a subsequence tending to a solution (A^*^, Aq) 
of the equation 

/-Xoa-f B(s,X{s,y))ds 
b{a,x,X{a,y))N'^{y)e " dyda. 

Defining N by p9|) . it remains to check that Aq > 0, which will imply J Ndxda < oo. Since ^ 
and A^*^ ^ 0, we have / N^{x)dx = 2 ff B{a,x)N{a,x) dx da > 0, so integrating equation d?]) in x, 
and in a from to A, one gets: 

A XM 

N{A.,x)dx + \Q j j N{a,x)dxda = J J B{a,x)N{a,x)dxda. 



Using the preceding calculation done for f No^^i^, x)dx and assumption ([^T|) . we deduce: 



A xm ^ / \ 

I f ff f ~\oA-J B[s,X(s,x))ds \ 

Xq N{a,x)dxda= / lB{a,y) — 2 / b{a,x,y)e " dx \ N{a,y)dy da > 0, 





which implies Aq > 0. To prove uniqueness, and to find a solution for problem ([8]), the proof is identical 
to the one of part 12.31 so we need to prove / Ncj) dx da > 0. 

Since N ^ 0, there exists yi €.]Q,xm[ such that A^(0,yi) > 0. Formula (fT9]l and assumption ([l3l) 
implies then N(0,x) > for x ^ X{Ai,y). Recursively, it stands for x ^ X'^(Ai,y) for all n € N. 
Under assumption (I42|) . X is increasing, so the sequence X"{Ai,y) tends to a limit / that verifies 
X{Ai,l) = I, so I = Xm and A^(0 < x < xm) > 0. Hence, f N(j)dxda > 0. 

3.2 Weak Theory When b{a, x, y) is Continuous in the x— Variable 

In this section, we extend the previous results to a larger class of parameters. First, we can relax 
assumption ([27|) . and replace it by the two following conditions. 

Coix)= [ U{a,x,X{a,y))e_Jl^^'^}^f^dady(^C\[{),XM]). (45) 

J " — T~^, 

birth term *erm 
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This assumption is necessary to obtain compactness, independently of e ^ 0, and to apply Ascoli 
theorem in E = C^(^[0,xm])- 

poo 

/ B{s,X{s,y))ds>\n{2), Vye]0,XM[. (46) 
Jo 

This condition means that there is enough birth along the characteristic curves. We check as in 
paragraph 12.11 that it is verified by the given examples of part II. 4i It is used to prove, with the 
preceding notations, that the eigenvalue //o,e > 1 (as in paragraph 13.11 it is no more necessarily equal 
to 2) and hence that the eigenvalue Aq > 0. 
We also relax the regularity assumptions on the coefficients, and we suppose: 

s G L^,(M^), re W^^°^{W\). (47) 

We obtain the following theorems. 

Theorem 3.4 Under assumptions (C^, (f-<^5D -( f77[ ), with the positivity assumptions 129\) and ^30\). 
there exists unique solutions (Aq, N, cj)) to the eigenproblems ^ and Xq > 0. Moreover, for all 
< ?7 < 1, Ne^oT" and (/)e^o^" G Li(M+ x [0,xm])- 

Theorem 3.5 Under assumptions U6\) . |^0[ ) and ( [^ ([^5|)-(f77|), with the positivity assumptions 
(4^ ) and [4^ , there exists unique solutions (Aq, N, (p) to the eigenproblems and Aq > 0. 

Moreover, for all 0< r] < I, Ne^o"^" G L^{R+ x [0,xm]) and (j) G L°°(M^), Suppa(t) C [0,^] x [0,XAf]). 

Proof. Refer to proofs of theorems EU [221 and EH 



3.3 Case of Equal Repartition After Division 

In this section, we consider the case 

h{a, X, y) = S^^RB{a, y), (48) 

With B G C(M^) and T verifying the Cauchy-Lipschitz conditions, B{a,x ^ xm) = 0. Our proof can 
be extended to the cases of unequal mitosis, i.e. if 6 is a sum of Dirac measures. We cannot use lemma 
11.11 since its formulae are not verified anymore, but the method remains. To make the details easier 
to understand, we also make the following assumptions (which is reasonable if we follow [6]): 

V0<a;^2f, Va^O, r(a,x)>0, 

(49) 

{{a,x) G M_|_x]0,X7v/[; s.t. r{a,x) = O} := xo{a) is an increasing curve. 

We make assumption (|lUp with T regular. It allows us to divide F by x and obtain a smooth function. 
We first establish the following lemma, which replaces lemma [TT] and on which is based the proof. 

Lemma 3.6 Let B G Cfc(M^), b defined by (JS^ and T G C-'^(M^) verifying assumption [J^ . If N is 
solution of problem ([7|) |j.^[ j then the two following identities stand: 

r -\oa-f{d^r+B) ( s,X{s,Y{a,2x))) ds 

A^(0,x) = 4 / B{a,2x)N{0,Y{a,2x))e o 'da, (50) 

a(x) 
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where a{x) is defined by a{x) = Inf[a ^ 0, Y{a,2x) ^ 

Mm(^,2x) f(a,2x) 



where f{.,2x) is the inverse function ofY{.,2x) and is so defined by f{Y{a, 2x),2x) = a. 

Proof. By the method of characteristics previously used, we obtain by straightforward calculation: 

r -XQa-f{dxT+B)(s,X{s,Y{a,2x)))ds 

N{0,x)=4 B{a,2x)N{0,Y{a,2x))e " da. (52) 



The flow Y(a,x) can be defined also by the following differential system equivalent to (fT2|): 



-fdxr(s,X{s,Y(a,x)))ds 

:^Y{a,x) = -T{a,x)e , a > 0, x ^ 0, ^^^^ 

Y{0,x) = x, x^O. 



For X ^ it implies that ■^Y{a, 2x) ^ so Y{a, 2x) ^ 2x ^ xm- Since we are looking for solutions 
verifying ()14p . formula (jSOh implies that A^(0,x ^ ^) = 0. Hence, in formula (j52p the integral can 
be reduced to a such that Y{a^ 2x) ^ and we find formula (f50]) . 

For X ^ 2|£^ assumption implies that ^(a, 2x) is decreasing with the a variable, so Y{a, 2x) ^ 
a(x) = 0. We make the change of variables a ^ a' = Y(a, 2x). When a = 0, a' = 2x and when a = +oo, 
a' tends to a limit I ^ 0, since for / > one has T{a,l) > so ^Y{a,l) < 0. Formula §3) comes, 

- J da:T[s,X{s,Yia,2x)))ds 

since we have da = ^y(a, 2x)(ia' = — r(a, 2x)e " da'. 

For ^ < X ^ since 2x > for a smah we have Y{a,2x) ^ so a(x) > 0, and the 
definition of a{x) implies r(a(x),2a;) ^ 0. Because we made assumption (jM]), T{a > a{x),2x) > so 
Y{a > a(x), 2x) is decreasing and we can define once more the change of variables a ^ a' = y(a, 2x). 
When a = a(x), a' = ^ and when a = +oo, a' tends to a limit / ^ 0. Formula (I5ip then comes. 

We now formulate an equivalent condition of assumption ([271) . used to obtain compactness: 

I' -/#-(r+B) s,X(s,y(a,a;)) 

h{x):= J e ^ V(a,.)G[o,^]^«eC^([0,XM]). (54) 

a{x) 



Formulated in terms of partial differential equations, it means that the solution v to (|28p satisfies 
V G C^([0, XAf], L^(IR+)) . For the adjoint problem, we replace assumption (fT6]) by: 

B{a,x^ X{a,XM)) =0. (55) 

To have uniqueness, we make the following assumption: 

VxG]0,xm[, VaGM+, B{a,x)>0. (56) 
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Theorem 3.7 Under the assumptions of lemma \3.6\ and the supplementary assumption (5^, there 
exists a solution {N,X) to problem (f7[) ((j^[ j and Xq > 0. Under the supplementary assumption (155)). 
which implies [T4\ ), there is a solution (p to the adjoint problem ^ related to an eigenvalue Xq > 0. 
Under the positivity assumption i56\). {Xo,N,(f)) are unique. 

Proof. We follow the same way than for the proof of theorem 12.11 for instance, so we detail only the 
specific points. The main new difficulty is that on formula (j5ip . we have divided by F (/(a, 2a;), 2a;) 
on a domain where it is strictly positive, so its inverse is well-defined, but it vanishes if x tends to 
or It does not matter however, since when it vanishes, its inverse is multiplied by N{0,a) which 
also vanishes. We are then led to use either of formulae (j50p or (j5ip according to where x stands. 



We first define a regularised operator Q^^\ : X ^ X on X = C{[0, ^]) by either of the equivalent 
formulae: 

OO / \ 2 

r -\-f{da:r+B)[s,X(s,Y(a,2x)))ds f 

Qe,\{9){x)=^ / B{a,2x)g{Y{a,2x))e « ' da + e g{x)dx 

a{x) 

Min(^,2x) f(a,2^) ^ 

B(f(a,2x),2x) ^ -\f{a,2x)- J Bis,X(s,a))ds I' , , 

r(/(a,2x),2x)^ ^"^^ ° da + ey,(.)dx. 



The difference with the previous regularisations of parts [2] and 13. II stands in the fact that we regularize 
differently: changing B in B^ would not change the value of Q{g){x = 0) which would remain zero. 

Lemma 3.8 Under the assumptions of lemma \3.(^ and the supplementary assumption ((5^[ ), for all 
e ^ 0, A ^ 0, the operator \ : X ^ X is compact. 

Proof. It suffices to take e = A = 0. The assumptions ensure that the operator Q = Qq^q is 
well-defined. To apply Ascoli theorem, let r/ > arbitrarily small, we look for > such that 

V5 G X, ll^lloo ^ 1, Va;i, xa G [0,^], < X1-X2I < i/^ A(xi,X2) = \G{g){xi)-G{g){x2)\ < V- 

We distinguish three cases: around around 0, and on the compact subset [5, ^ — 5] with 6 > 
small enough. 

1. For xi, X2 close to we use formula (150p to define Q, and remark that lim^.^£M a{x) = 00. 
Indeed, for x close to F(0,a;) < and under assumption ()49p the function F(a,x) remains 
negative till a = ao{2x), ao{2x) defined by F(ao(2x), 2x) = 0. The curve Y{a, 2x) defined by ([53P 
increases for a ^ ao(2x), so for a ^ ao(2x) one has Y{a,2x) ^ 2x > ^ : hence, a{x) > ao{2x). 
Under assumption ([39]) it is clear that limax-^^^^^ ao{2x) = 00, so lima^^xM Oi{x) = 00. 

Assumption (|54l) then implies lima^^xM ^(/)(^) = uniformally in g, if \\g\\oo ^ 1- 

2. For xi, X2 > close to 0, we have G{g){0) = since B{a,0) = 0. Since B is uniformally 
continuous, and the integral operator is uniformally convergent thanks to assumption (j54p . we 
can bound G{g){xi^2) uniformally for xi^2 small enough and H^Hoo ^ 1- 
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3. For 5 < xi, X2 < ^ — S, we use formula (jSip to define Q. This formula gives a classical form of 
Q as an integral operator: we know that it is compact as soon as the kernel under the integral 
is continuous and bounded. It remains to prove that T{f {a, 2x),2x) does not vanish. 

For X < ^^sM^ has 2x < xo(0) ^ xo{a) for all a ^ 0, so under assumption (jl9]) we have 
r(a, 2x) > for all a > 0. 

For ^ ^ X < 2f , for all a we have proved above that /(a, 2x) ^ a{x) > ao{2x), which implies 
by definition of ao{x) that T{f {a, 2x),2x) > 0. 

Since r(f{a,2x),2x) > on the compact subset {x,a) G [|, - |] x [0,Min{2x, ^)], and F 
is continuous, it reaches its minimum F*"-^ > : this ends the proof of lemma 13. 81 

Lemma 3.9 For all e > and A ^ 0, under the assumptions of lemma \3.8l there exists a unique 
l^e,\ > o.f^d a unique ^ E X such that Ge,\{^£x) = /^e,A-^£ A' ll^e aII^ ~ Moreover, X ix^^x 
is a continuous decreasing function, with fie^oo = and Hefi = 2 + e^. 

We let the reader check the proof (equivalent to that of lemma [23]) . We define, as soon as < 1, a 
unique > such that fj^c^x^ = 1. We denote as before the associated eigenvector with HA'^PHx = 1- 
As in lemma [231 it comes from the previous study that the family (A'^)o^e^i is compact, so we extract 
a subsequence tending to (Ao,iV'') G M_|_ x E. Under assumption ([56j) . formula ([5T|) implies that if 
N^{xi) > then iV°(x g]^, 2f [) > 0, so recursively it implies 7V°(]0, > 0. 

The resolution of the adjoint problem is made as in theorem l2.2[ The solution of the adjoint problem 
([8]) can be written now as 

oo a , V 

f X(n r) -Xoa- f B{s,X{s,x))ds 

0(O,x) =2 / i?(a,X(a,x))0(^^^^)e i ^ ' da. 



Since 7^ and ^ there exists < xi < ^ where (p{xi) > 0. So // Ncpdxda > 0, which 

implies Aq > 0. The proof of uniqueness of A'^ and cp is the same than for theorem 12.31 



4 Extensions 

4.1 Resolution of a Model With Multiple Cyclins 

As already mentioned, there is a whole variety of proteins and cyclin/CDK complexes which play a 
role in the cell cycle, and we can also want to structure the model by the DNA content or the size of 
the cells, etc. Hence, that would be useful to include in the model the action of several variables. 
Let us suppose that we have n variables playing a role in the cell cycle, and denote them by Xi with 

/ n 

1 ^ i ^ n. We write x = (xi, F = (Fi, F„), \x\ = J^xf, 

and we define an order on M" thanks to the cone by x ^ y <^=^ Vl^i^n, Xj^yj. 
The model ([7]) may be generalised by: 



' -§^N{a, x) + Vx • [F(a, x)A^(a, x] + (Aq + B{a, x)) N{a, x) = 0, a ^ 0, x ^ 0, 

Af(a = 0,x) = 2/6(o,x,y)A^(a,y)(ia(iy, (57) 
N^O, fNdadx = l. 
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In the same way, the adjoint problem becomes: 



-§-J - r(a, x) • V^<t, + (Ao + B{a, x)) <P = 2 J 4>{0, y)b{a, y, x)dy, a ^ 0, x ^ 0, 

(58) 

' ^ 0, / (l)Ndady. = 1. 
We use the notation V = (rj)i^j^„ and suppose that rj(o, x) = XiTi{a^K), with 

f(a,x) >0 forxPsO, 

(59) 

r(a,x) ^0 for X ^ xa/ = {xim, ■■■,XnM) 

Concerning division, we have the same relations ^ and ([3]) than previously seen. Characteristics and 
their inverse flow are still defined by (jl2p and (|13p . Assumption (j59p implies that for all ^ x ^ xm, 
X(a,x) ^ xjv/, so we add the same conditions (fH|) and ([TUD to problem ([57]) . 

n 

Theorem 4.1 Replacing x, y £ M-|- fey ^ ^'^^ Y ^ 1^+; replacing ^ 6?/ V • F, anrf y'^'i 6?; Y\ vT 

^ i=l 
7j ^ 1, a// t/ie preceding theorems \2.1\ \2.2l \2.3l \3.1\ and \3.7\ extend to the n— dimensional case under 

the equivalent assumptions in n dimensions. 

4.2 Asymptotic Behaviour of the Linear Evolution Problem 

Having solved the eigenvalue problem, we are now able to characterize the asymptotic behaviour of 
the solution n to the time-dependent problem ([T]). First, we establish a General Relative Entropy 
Inequality, using the same formalism than in [27] and ^28j . 

Proposition 4.2 Let p(t,a,x) ^ 0, n(t,a,x) smooth solutions of problem ([T]) and <I>(t, a,x) ^ 
smooth solution of the adjoint problem: 

- - r(a, x)-§^<^ + B{a, x) $ = 2 / $(t, 0, y) b{a, y, x) dy, a^O, x^ 0, 

(60) 

$ ^ 0. 

Then we have, for any function H : 



^(f )(0 = i I Pit, «' ^Mt, a, x)tp^t^a,.)^oH{^^^)dadx = 2 /p(t, a, x^t, 0, y)b{a, y, x) 

dadxdy. 



TT( n{tfl,y) \-r. TT( n{t,a,x) \-,, , , i ij/ / w(t,0,i/) \ / n(t,a,x) n{t,0,y) 



i) If p > never vanishes, for H convex, we get ^T~(-{^){t) ^ 0, 

and if H is strictly convex, 4[Hi-){t) = iff n = Cp on Supp{b) with C ^ constant. 



dt ^p ' 

ii) If p ^ 0, for H convex, positive and non- decreasing, we also get ^'H{^){t) ^ 0. 
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Proof. If p > everywhere, Ip^o = 1 ^^id we have by straightforward computation: 

ft T) (/ 71 ff T) 71 I 

^(p^H{-)) + {p<^H{-)) + [T{a,x)p^H{-)] = -2pH{-) / ^{t,Q,y)h{a,y,x)dy. 
ot p oa p ox p p J 

Integrating this identity and denoting 'H{^){t) = J p{t, a, x)^{t, a, x)H {^^^^^^^dadx, we get: 
- = 2 y a, x)<I>(t, 0, y)Ka, y, x) M:^) - ^ (^) 

because we have noticed that 

p{t,0,y)' ^p{t,a,x) p{t,0,y)^ 

Up^O, we replace by H{^)lp^o in the calculation, since 9;,(p(/>ii'(^)]lp^o) = Hp^o^^P^^^lf )) 

and 

y ^ (r^^ [rm:^ " M^^'^^ ' ^ ' ' ''^ ^' " ^ 



Theorem 4.3 Zet no € L^i^'j^(f){a,x)dxda), Supp{nQ) c M4. x [0,xm]) T satisfying the Cauchy- 
Lipschitz conditions and B £ L°°(M^). 

There exists a unique solution n € C(M+; L"'^(M^)) of the following problem: 

' §-^n + -§^n + -§^[T{a,x)n] + {B{a,x) + Xo)n = 0, a^O,x^O, 
< n{t,a = 0,x) = 2 J b{a,x,y)n(t,a,y)dyda, (61) 
n{t = 0,a,x) = n^{a,x). 

Moreover, we have the following inequalities, if\n^{a,x)\ ^ CQN{a, x), 

(i)yt^O, \n{t,a,x)\ f^CoN{a,x), 

(a) ^ n2 =^ ni{t,x) ^n2{t,x), 
00 00 

(Hi) J n{t,x)4>{a,x)dxda = J n^{a,x)4>{a,x)dxda, 



00 00 

(iv) J \n{t,x)\4>{a,x)dadx = J \n^{a,x)\4>{a,x)dxda. 


Theorem 4.4 Under the assumptions of theorem \4-3[ and if we suppose also that satisfies 

!n°(a,x)| ^iV(a,x), \dan^ia,x) + d^{rn^){a,x)\ ^CiiV(a,x), 

the solution to [61\) satisfies also \dtn{t,a,x)\ ^ (Ci + Aq + ||-B||L°°)-^(a,a^)- 
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Theorem 4.5 Under the assumptions of theorem \4 -SI and either those of theorem \3.4\ or of theorem 
\3.5l defining mP = J n^{a,x)4>{a,x)dx da, the solution to ( fg7]) satisfies: 

I \nit,a,x)-m<^Nia,xMia,x)dxda 10 as t ^ oo. 

Proof. Since the proofs are the adaptation to our model of those of theorems 3.1, 3.2 and 3.4. of 
|31j or of theorems 4.3 and 3.2 of [28j, we let the reader check them. 



4.3 Application to a Two Phase Model 

As previously seen in part 11.11 equation ([7]) can be considered as a simplification of the linearised 
eigenvalue problem for the two-compartment model described in [6]: 



^P + -^[T{a,x)P] + {\ + B{a,x) + di + L{a,x))P-GQ = 0, a^O,x^O, 
{X + G + d2)Q = L{a,x)P, a^O,x^O, 

(62) 

P{a = 0, x) = 2 J b{a, x, y)P(a, y)dy da, 

P, Q ^ 0, f P + Qdxda = 1. 
The adjoint problem is: 

' -^</. - T{a, x)^ + {X + B{a, x) + di + L{a, x)) 4> - L{a, x)i; = 2 J (f>{0, y)b{a, y, x) dy, 

{X + G + d2)^ = G(p, a^O, x^O, (63) 

, (/), V ^ 0, f(t>P + i;Q dx da = 1. 

Here P and Q denote respectively the proliferative and quiescent populations of cells, G is the recruite- 
ment function, L(a, x) the number of cells going from the proliferative to the quiescent compartment, 
and di, d2 are the death rates of each population. For the sake of simplicity, we have limited our study 
here to the case when L is constant, and we make the following assumptions for the coefficients. 

di>0, d2^0, L{a,x) = L^O, G ^ 0. (64) 



Theorem 4.6 Under the assumptions of theorem \2.S\ and with the supplementary assumption ^64^ , 
there exists unique solutions {P,Q) € E'^, (0, -0) G E'^ to problem and l63\) for a unique A € M. 
Moreover, denoting Xq the eigenvalue of problem ^ and we have: 

Ao = A + di + L ^f^^ > 0, (65) 
X + G + d2 

or also, defining G+ = G + d2, d+ = di — Aq and L+ = L + di — Xq 



^ ^ -{G+ + + V{G+ + - 4((i+G+ + d2L) ^^^^ 
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The following estimate stands for A, with Li = L + di : 

A > A := ^ + Li) + + Li)2 - 4(diG+ + Lds)^ . (67) 

Proof. We reduce the system to a single equation on P, since we can write Q = ^^q^ P- Writing 
Ao = A + di + L^±^, we find: 

^^ + ^[r(a,x)P] + (Ao + 5(a,a;))P = 0, a^O, x^O, 

P(a = 0,x)=2j b{a, x, y)P{a, y)dy da, (68) 

^^0, jp{i + ^)dxda = l. 

We only have to change the normalization and the previous study applies to this case: theorems 12.11 
to 12.31 give us a unique solution (P, Q, cj), ip, Aq > 0). We denote the eigenvalue Aq = /(A) where / is a 
continuous increasing function with a unique singularity for A = — G+. Since Q = x+g+ need to 

have A > — Ai, so for each Aq G M, there exists a unique convenient A g] — G+, +oo[ which is given by 
formula (j66p . The term under the square root is always nonnegative since 



(G+ + L+f - Ad+G+ + d2L) ^ (G+ + L+f - 4.L+G+ = (G+ - L+f ^ 0. 
The inequality (136]) is given by the fact A > A2 where /(A2) = 0. 

Discussion. Prom theorem 14. 6l we deduce as in theorem 14.51 the long-time convergence of the solution 
of the linearised problem ([6]) towards {P,Q)e^^. However, what is experimentally observed is either 
convergence towards a steady state or exponential growth but only in the early stages (Gompertzian 
growth: cf. [20] and references therein). In [8] polynomial infinite growth for 15 cell lines is shown, and 
in [17] a single-cell model is built, which is able to exhibit such behaviours. The linearised problem 
cannot take into account such phenomena, due to feedback answer or saturation effect: it can only 
come from a non-linearity of the model. 

But if the linear renewal equation has a relatively simple asymptotic behaviour, the theory for 
nonlinear models is much more complicated. Several behaviours are possible: chaotic, periodic, con- 
vergence towards stable steady states (for recent references on nonlinear population models, see for 
instance [3], 0, [TO], [H], [TO], [TO], [25], [29] or [33]). 

In [6], following [20] and |21j . it is proposed that the non linearity comes from the term G{N{t)) where 

00 00 

the weighted total population A^(t) is defined by A^(t) = / J[(t)*{a,x)p{t,a,x)+il)*{a,x)q{t,a,x)]dxda, 



where cj)* , '0* are given weights. The recruitment function G is taken equal to 

To study the behaviour of the model, the method of [6] is inspired of the principles of General Relative 
Entropy. It is based on estimates and on the construction of subsolutions and supersolutions of a 
quantity 

00 00 

S(t) = j j [(j){a,x)p{t,a,x) + 'il){a,x)q{t,a,x)\dadx, 
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where {(p, ip) is the solution of the adjoint linearised eigenproblem (j63p for a proper value of G = G{N) 
(see also |12j for the application of this method to another model). Proposition 2.7. of p] shows that 
unlimited growth can be obtained under the two following conditions. 

(H9) For all N < 00, the eigenvalue A(A^) corresponding to G = G{N) in ([62|) satisfies A(A^) > 0. 
(HIO) For each corresponding solutions to the systems (f63]) with G = G{N), denoted {(p^jipiy), there 
exists a uniform constant C„ such that cp* ^ Cu4>n and ^p* ^ GuipN- 

Proposition 2.5. of [6j shows that subpolynomial growth can be obtained under the following con- 
ditions, if (i2 > and 02 > 0. 

(H7) For G = G{oo) = 02 > 0, the first eigenvalue of ([62]) is A(cx)) = 0. 

(H8) For the corresponding solutions to (f62]) and (f63l) . denoted respectively {P2,Q2) and (02, "02)1 
there exists positive constants C2 and C3 such that C3(/)2 ^ 0* ^ C'2(/>2 and (73-02 ^ "0* ^ G2ip2- 

But to obtain exactly (H7), it has been observed in the numerical simulations of [6j that all the 
parameters di, d2, L, 02, ai had to be related and chosen very carefully: a very small change in one 
of the parameters implies that X{oo) 7^ 0, and the system either is bounded (if A(oo) < 0) or grows 
exponentially (if A(oo) > 0.) This can also be seen by taking a closer look to formula ()66p : indeed, we 
can also write it as 

A = -KG.d. + Ld2) 

G++L+ + V(G+ + L+)2 - AG+d+ 
Since the denominator of this formula is always positive and bounded, one has X = ijf 

G+d+ = -Ld2. (71) 

The simulations carried out in [6j were all done with ^2 > 0, so this formula is verified punctually, for 
special values of the coefficients G, d2, di, ^2 and L linked by ([7T]). 

From a biological point of view, this obligation to have coefficients linked by such a relation seems 
hardly justified. But we can also assume ^2 = : from a biological point of view it can be verified for 
some kinds of cell populations - stem cells for instance: in the quiescent compartment indeed, there is 
no reason why the cells should die (see [T] or [20], and the references therein). 

In this case, and supposing also that lim G{N) = (which is indeed realistic) we see that condition 

(|7ip will be always verified: we can now obtain a "robust" subpolynomial growth - I call "robust" 
a subpolynomial growth which remains true for a whole range of parameters di, L, n, ai. This is 
expressed by the following proposition. 

Proposition 4.7 Let us suppose d2 = 0, G{N) defined by i69\) with 0,2 = 0, T{a,x), B{a,x) given 
functions verifying the assumptions of one of the theorems \2.3\, \3.1\ or \3. 7[ We denote Aq > the 



£P2 + -^[T{a,x)P2] + {B{a,x) + di + L) P2 - Q2 = 0, a ^ 0, x ^ 0, 



eigenvalue of p^, and suppose that < di < Aq and that L > \q — di. The case lim G{N) = G = 

N^oo 

can be represented by the following system: 

da 2 ~r 

Q2 = {L + di- Xo)P2, a^O, x^O, 
^2(0 = 0,x) = 2/ b{a,x,y)P2{a,y)dy da, 
I P2, Q2>0, JP2 + Q2 dxda = 1. 



(72) 
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The adjoint problem is: 



-^(/)2 - r(a, x)^+ {B{a, x)+di + L) </.2 - L^s = 2 / (/>2(0, y)h[a, y, x) dy, 
Lij2 = {L + di- Ao)(^2, a ^ 0, X ^ 0, (73) 

02, V'2 ^ 0, / (/>2P2 + V'2Q2 C^xda = 1. 

Under assumption (HIO), and under assumption (H8) adapted for the prohlems(72) {73), there exists 
a constant C > such that: 

N{t)i^Ct^, lim iV(t) = +00. 

t—l-OO 

Proof. Since di < Aq formula (j70p with d2 = implies that (H9) is verified, so proposition 2.7 of [6] 
can be applied and proves unlimited growth. 

The proof of the subpolynomial growth is based on the same tools than in [6], proposition 2.5: the 
only difference is that since d2 = 0, at infinity we have A(cxd) = G{oo) = 0, so the equation for Q 
in (I62p only expresses that Q tends to infinity whereas (A + G)Q remains finite. We will obtain a 

relevant problem by an asymptotic analysis: formula ()70p can be written, if G, A ^ 0, A ~ L+di -\l ■ 
We can replace it in the second equation of (f62|l and find: GQ i^J'i-Xo ~ ^otiiig Q2 = GQ, we 
obtain a problem for the couple (P2, Q2) which remains meaningfull if A, G vanishes, and by choosing 
an appropriate normalisation its limit is (1721) . which adjoint is (j73p . 



We define S2{t) = j J [(j)2{ci', x)p{t, a, x) + ip2{a, x)q{t, a, x)]dadx, and we calculate 


dS2{t)^^^ = [[L^[.(^L + B + di)p + G{N{t))q - ^ - ^(Fip)] + MLp - G{N{t))q) } dxda, 



dt J J { da dx 

dS2it) f f ( . , X , / u N , , f f 



(t) = G // {(l)2{a,x) — ^2ia,x)^q{t,a,x)dxda = G // - — — ^—ip2qdxda^CGS2{t). 



dt J J J J L + di — \q 

Since N{t) ^ C^S2{t) thanks to (H8) we obtain 

dS2{t) a,e- 

^ ^^^^')^" + (C352(t))" ^^'^ 

Since it has been proved in p] that = a(t + to)" 1 for a large enough, is a supersolution to ([H]), so 
for to such that S(0) ^ 5*2(0), we have by the comparison principle S2{t) ^ 5](t). It ends the proof. 

For the numerical simulations (see figure [2]) , we take the same values of the parameters than in [6] , 
except that d2 = : we define T{a,x) by ([22]) with ci = 0.1, C2 = 0.075, n = 3, C4 = 0.4, r2 = 1.95. 
We take ai = 8, 9 = ln=^ with A: = 1, 2, 3 in the definition (j69p . We define L{a, x) as in [6j by: 

L(o, x) = A3^^j—^\A,+oo[{a), with 72 = 5, A3 = 4 A2 = 2 A = 18. 

We define B{a,x) by S(a,x) = ^4r^l[^*,oo[(a) with h = 1.2 /c2 = 1.5 71 = 5 A* = 23. We 
obtain polynomial growth by taking for instance di = 0.01 and it remains true if we make small 
changes of any coefficient. But if di becomes too big, for instance if we take di = 0.05, we obtain 
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n=1, P{t)m(t) 
-- n=1/2, P(t)/N(l) 
n=1/3, P(t)/N(l) 
n=1,P(l) 
n=1/2, P(t) 
n=1/3, P(t) 





3000 4000 5000 6000 7000 8000 9000 



Figure 2: Evolution of the total cell population jj q{t,a, x) dx da + // p{t,a, x) dx da for a tumoral 
tissue. Left: with di = 0.01, with different values of n = 1 (lower solid line curve, left), n = 1/2 
(medium solid line, left), n = 1/3 (upper solid line, left) at a Log-Log scale, for N{t) (solid lines) P{t) 
(dotted lines) and P{t)/N{t) (dashed lines). Right: with di = 0.05, we see exponential decreasing 
(Log scale, the three curves for N{t) are superimposed). 



solutions exponentially vanishing. This is explained by proposition 14. 7t indeed, when di increases, it 
becomes bigger than Aq and formula ([70|) shows that it implies A(A^) < for all N. 



Conclusion. In this article, we first solved the eigenvalue problem ([7]) under fairly general assump- 
tions, provided the continuity of the repartition function b or, in the case of equal repartition after 
division, of the birth rate B (the generalisation to coefficients is a work in progress). Using General 
Relative Entropy Inequality, we deduced from it the asymptotic behaviour of problem ([T]). 

We then applied these results, in part 14.31 to the study of a non linear two cell-compartment model 
given by equation ([6]), model which was first introduced by F. Bekkal Brikci, J. Clairambault and B. 
Perthame in [6j and [7] to study the action of proteins on the cell cycle (it can also be considered as 
a generalisation of the pure size-structured two-compartment models studied by M. Gyllenberg and 
G.F. Webb in [20] and [H].) 

Finally, we exhibited a case of "robust" polynomial growth, which reveals coherent with the results 
of d] and [IT]. 



This last result could lead to two biological interpretations. First, when the population becomes 
larger, the formulation of the eigenvalue problem (f72l) (f73l) . where Q had to be replaced by Q2 the 
limit of G{N)Q, seems to show that the number of quiescent cells tends to infinity more rapidly than 
the number of proliferating cells: indeed, it seems that -P(i) ~ G{N{t))Q{t) ~ so the relative 

number of proliferating cells, given by R{t) = ^py, seems to vanish like It is coherent with the 
results of [20] and it is also confirmed by numerical tests (see figure El left side) . From a biological 
point of view, it seems true in many cases, for instance for stem cells (see [1]), that most of the cells 
are quiescent. Second, this qualitative study seems to emphasize the crucial importance of apoptosis 
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for homeostasis or tumour growth (see [23]): indeed, if di = Xq is reached, the number of tumour cehs 
would decrease rapidly instead of growing to infinity. 
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